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Motivated in part by quantum criticality in dissipative Kondo systems, we revisit the finite-size 
scaling of a classical Ising chain with l/r^~^ interactions. For | < e < 1, the scaling of the dynamical 
spin susceptibility is sensitive to the degree of "winding" of the interaction under periodic boundary 
conditions. Infinite winding yields the expected mean-field behavior, whereas without any winding 
the scaling is of an interacting cj/T form. The contrast with the behavior of the Bose- Fermi Kondo 
model suggests a breakdown of a mapping from the quantum model to a classical one due to the 
smearing of the Kondo spin fiips by the continuum limit taken in this mapping. 
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Quantum criticality describes the collective fluctua- 
tions of a continuous quantum phase transition. The 
standard theoretical picture [1] is that quantum critical 
fluctuations in d spatial dimensions can be described in 
terms of the fluctuations of an order parameter (a classi- 
cal variable) in d-\- z dimensions, where z is the dynamic 
exponent. However, there has been much recent interest 
in quantum critical points (QCPs) that do not conform 
to this quantum-to-classical mapping [2]. 

One setting to explore this issue is quantum impurity 
models that involve one or more dissipative bosonic baths 
having a sub-ohmic (e > 0) power-law spectrum 



^[(5(cc; - Wp) - 5{u; + Wp)] 



^ ^sgn(cj)6(cc;c - \u;\). 

(1) 

Here Wp is the dispersion of the bosonic bath(s) and O is 
the Heaviside function specifying the high-energy cut-off 
ujc- The quantum-to-classical mapping leads to a clas- 
sical Ising chain with l/r^~^ interactions, which has as 
its continuum limit a local (f)^ theory in + 1 dimension 
[3]. For 1/2 < e < 1, the 0^ theory has a Gaussian fixed 
point with a violation of uo/T scaling. However, the QCP 
in a large-A^, fixed-/^ limit of a spin-rotation-invariant 
SU(7V)xSU(/^7V) Bose-Fermi Kondo model (BFKM) sat- 
isfies co/T scaling and, hence, is interacting [4]. Nu- 
merical renormalization-group (NRG) studies have also 
found interacting QCPs for 1/2 < e < 1 in two quantum 
models with Ising anisotropy: the dissipative spin-boson 
model (H and the Ising BFKM 0. The failure of the 
quantum-to-classical mapping has been attributed in the 
spin-rotation-invariant case to the Berry phase 0, but 
remains to be clarified for Ising symmetry. 

This Letter reports a two-pronged study of the classi- 
cal mapping of the quantum-critical Ising BFKM. First, 
we show that, for 1/2 < e < 1, the scaling of the dy- 
namical spin susceptibility of the classical model as a 
function of system size and wave vector (corresponding 
in the quantum problem to temperature and Matsubara 
frequency, respectively) depends crucially on the "wind- 



ing" of the interaction to extend the range of the interac- 
tion beyond the finite system size. For infinite winding, 
one recovers the mean-field a;/T^/^^~^^^ scaling expected 
from the corresponding (j)^ theory, in which the quartic 
coupling acts as a dangerously irrelevant variable. Lim- 
iting the winding introduces an additional temperature 
dependence dominating that from the dangerously irrel- 
evant coupling and giving rise for zero winding to uo/T 
scaling. Second, we provide strengthened evidence for 
non-mean-field behavior in the BFKM. These findings 
can be reconciled through a breakdown of the quantum- 
to-classical mapping arising from failure of the continuum 
limit to fully preserve the effects of Kondo spin fiips. 

Bose-Fermi Kondo model and the classical Ising 
chain. — The BFKM describes a quantum spin coupled to 
both a fermionic band and a bosonic bath. The Hamil- 
tonian for the Ising-anisotropic case is 



-^P ^po- ^pa 



Wbfkm = Jk^ • Sc 

^P^l^p^ (2) 



where S is a spin-1/2 local moment, the c^cr's represent a 
fermionic band with a constant density of states 6{uj — 
Ep) = No and on-site spin Sc, and the (j)pS represent a 
bosonic bath having the spectrum in Eq. ([T]). For fixed 
NqJk, a QCP g = Qc separates Kondo (or delocalized) 
and local- moment (or localized) phases. For ^ = 0, the 
BFKM reduces to the conventional Kondo model. 

The partition function associated with Hbfkm is [1] 
^BFKM ^ Trexp(-6'(^p), where 

Xo^(r-rO-/Ce(r-rOl| (3) 



and the trace is over the spin degrees of freedom. Xo ^ 
t'\ with short-time cutoff Tc = 2Ti/ui)c^ encodes the bath 



2 



spectrum. /Cc(r — r'), which has a Fourier transform 
JCc{iuJn) = i^\^n\ with = ttNqJk^ comcs from inte- 
grating out the fermions. 

Trotter decomposition of «S/^p (with time sUce tq) re- 
quires identifying the effect of the transverse-field term 
r5'^ with that of a nearest-neighbor term in the corre- 
sponding Ising model described [9*, ^lO^ by 



Z ^ Ti exp 



L 

i=l 



(4) 



The nearest-neighbor interaction i^nn = — 21n(ror/2), 
while the long-range interaction K\^{i — j) results from 
discretizing Xo^i'^ ~ '^') — ^c{^ — t'). The bosonic spec- 
trum specified by Eq. ([1]) gives rise for max(ro, Tc) <C r < 
j3/2 to a simple power-law behavior along the imaginary- 
time axis: 



Xo'(t)-1/M = 



for < e < 1. 



(5) 



The inverse temperature (5 = 1/T in the BFKM sets 
the length L = I3/tq of the Ising chain, with the pe- 
riodic boundary condition 5'£^2 — enforced by the 
trace operation. Varying the bosonic coupling g amounts 
to changing the effective temperature T of the classical 
model (which bears no relation to T in the BFKM). 

Classical Ising chains with this type of long-ranged 
interaction have been studied for decades, and it is 
well established that the phase transition for e = is 
Kosterlitz-Thouless-like [ol, [llj and is described over the 
range < e < 1 by a local (j)^ theory [1, [l2[. For 
< e < 1/2 the phase transition is controlled by the 
interacting Ginzburg- Wilson-Fisher fixed point, whereas 
for 1/2 < e < 1 mean- field behavior obtains. 

In general, if the classical system is below its upper 
critical dimension, then close to the critical tempera- 
ture Tc the only relevant scale is the finite system size, 
and the static part of the order-parameter correlation 
function should scale as Xstatic(^) = ^^^(^/^oo), where 
^oo Static (^ = oo) ~ (T - %y . Combined with 
Xstatic(^ oo) (T -Tc)~^ , we have x ^jv. From 
the hyperscaling relation 7 = — r\)v and the fact Q 
that 77 = 1 — e, we end up with x = 1 — e: at T = 7^, 



Xstatic(^) ~ 



for < e < 1/2. 



(6) 



Above the upper critical dimension, the finite-size scaling 
is complicated by the presence of dangerously irrelevant 
variables that introduce additional scales to the problem 
and destroy hyperscaling [13]. Extensive theoretical work 



(7) 



la hj has concluded that at T = %, 



Xstati 



,{L)r^L'l^ for l/2<e<l. 



Beyond the static limit, the full scaling properties of 
the susceptibility have not been clarified. We are aware 
only of a finding (Fig. 4 of [l5|) that x('?", L) fails to show 



any scaling collapse in terms oir/L. Separately, a study 
of the phase diagram has hinted at the importance of 
the winding of the interaction around the finite-size Ising 
ring [1^. These considerations motivate us to carry out 
a careful analysis of the scaling of x with particular at- 
tention to the effects of winding. 

Finite-size scaling of the classical Ising chain. — The in- 
teraction Xo^{^~^') is specified by Xo "^(^) = I^Xo "^(^ + 
iO+) = 7r^^cc;-^~^sgn(cc;), which follows from Eq. ([T]) with 
the cutoff ujc taken to infinity: 



Xo 



\t-t')= I 
Jo 



cosh [u{f3/2- |r-r^|)] 



duJXo (^) 



sinh(cj/3/2) 
Xo-'(r-r^T,oo), (8) 



where 

Xo'{T,T,M) = 7rr{2-e)g' J2 [{\r\ + np)' 



M-1 



■(2-e) 



n=0 



+ ((n+l)/?-|r|) 



-(2-e) 



(9) 



We study the Ising chain with effective spin-spin in- 
teraction Xo ^('^ ~ ^7 ^) fo^ different values of M, the 
maximum number of times that the interaction wraps 
around the Ising ring of length L. The \r — r'\ 00 
behavior is not altered by choosing a finite value of M. 
We use a cluster-updating Monte Carlo scheme [13, Q 
to measure the dynamical susceptibility xi^n^T) at cou- 
pling g = gc (i-e., T = 7^ in the classical model); 
cOn = 27rn/(roL), n = 0, ±1, ±2, . . . are the Matsubara 
frequencies of the quantum model. A binning analysis of 
our data shows that the relative error Ax/x ^ 
the integrated autocorrelation time is small. 

For < e < 1/2, the scaling properties do not depend 
on the choice of M and the static susceptibility obeys Eq. 
(j6|). As a result, the dynamic susceptibility scales as 



x(c^n,T)=T-(i-^)^(c^,/T). 



(10) 



Figure [H illustrates this uJn/T scaling for e = 0.4 with 
M = 1 [part (a)] and M = 10^ [part (b)]. 

For l/2<e<l, by contrast, the scaling properties 
depend crucially on M, as exemplified in Fig. [1] for e = 
0.8. The scaling plot for M = 1 [part (c)] implies that 

Xi^n) 



-0.2 



and x(^) ^ T-^ ' Only for M 



00 



do we recover the anticipated mean-field result of Eq. 
([7j), as shown in Fig. H^d) for M = 10^. Between these 
extremes, there is a slow crossover with increasing M, 
exemplified by the evolution of Xstatic(^) Difficulty 
in recovering mean-field physics in this model has been 
observed in the context of the phase diagram [l^, but 
the origin of this behavior and its implication for cOn/T 
scaling have not previously been recognized. 

In obtaining Fig. (TJ we set = and Knn = because 
the construction of clusters via continuous bond proba- 
bilities is efficient only for pure power-law interactions 
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FIG. 1: (Color online) Finite-size scaling of the susceptibility 
below (e = 0.4, upper panels) and above (e = 0.8, lower pan- 
els) the upper critical dimension, for tq = 1/16, n = Knn = 0, 
and winding numbers (a,c) M — 1 and (b,d) M — 10^. For 
e > 1/2, the temperature scaling depends on the choice of M. 
Similar results have been obtained for various tq between 1/4 
and 1/64. 





(a) 


: (b) 


10 


-- /(cc) = 0.85a;-°-2 


. "^^^ — fix) = 0.43cc-o-4» 


1 






10-1 










: . x(^n = 0,T) \ 


10-2 


. . xi^n,T = 1.22 X 10-^) \ 


- . x(u;^,r = 4.88 X 10-4)\ 
; - fix) = 0.81a;-o-2i \ 


10-3 






10-4 10-2 1 
or T/D 


10-2 1 
LOn/D or T/D 



FIG. 2: (Color online) Static and dynamical susceptibilities 
for e = 0.8 with k = 2 and Knn specified by F = 0.75 and 
TO 1/8: (a) M = 1, and (b) M = 5 x lO'^. 



17|. For nonzero and i^nn, the calculations can be car- 
ried out only for a more limited range of uOn and T, but 
the same behavior is found, as illustrated in Fig. O 

The difference between the cases of infinite and zero 
winding can be understood as follows. The infinite- 
winding interaction Xo'^i^n^T^ M = 00) depends only 
on cOn- The temperature dependence of xi^n^T) is then 
entirely determined by that of the "spin self-energy" 



7W(cu„,T) = Xo i^n,T) - l/xK,T). 



(11) 



The temperature dependence of M. is controlled by the 
dangerously irrelevant quartic coupling of the local (j)"^ 
theory, and is proportional to T^/^. On the other hand. 
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FIG. 3: (Color online) Static susceptibility at the QCP of the 
Ising BFKM for e = 0.8, calculated using the NRG with dis- 
cretization A = 9 and various values of the bosonic truncation 
parameter 8 < A/'^ < 20. See for further details. 



the zero- winding interaction XQ^ij^n = 0,T, M = 1) ~ 
T^~^, so x~^(^7^n) acquires temperature dependence 
from both Xo^^ and Ai. The leading term varies as T^~^, 
which overpowers the T^/^ from the dangerously irrele- 
vant coupling, so the uo/T form of Eq. (p!Q|) ensues. We 
have confirmed that this distinction between M = 1 and 
M = 00 is robust over the range 1/64 < tq < 1/4. 

Quantum critical behavior of the BFKM. — We have 
systematically extended previous NRG calculations 
to investigate the robustness of the quantum critical be- 
havior. Figure [3] shows data for e = 0.8 that can be fitted 
to Xstatic(^) ^ over more than 20 decades of T with 
X = 1 — e to better than 1% accuracy. The value of the 
exponent shows no systematic evolution as the bosonic 
truncation parameter Nb increases from 8 to 20. 

These NRG results, which are consistent with uo/T 
scaling, are in stark contrast with the cj/T^/^^"^^^ scaling 
and exponent x = 1/2 found in the regime 1/2 < e < 1 
of the classical Ising chain (with infinite winding). This 
disparity signals the failure of the quantum-to-classical 
mapping for the QCP of the Ising BFKM. 

The most likely source of the breakdown lies in the 
replacement of the spin-flip (transverse Kondo-coupling) 
term in the quantum model by a (9r^)^ term in the local 
0^ theory. In the BFKM, spin flips are essential for the 
formation of the Kondo-screened state and, by extension, 
central to the nature of the QCP displaying a critical 
Kondo effect. (Spin flips are the one-dimensional analog 
of the vortices in, for instance, the short-range XY model 
in two dimensions.) In the local (j)^ theory, however, the 
(drcj))'^ term (i.e., the oo'^cj)'^ coupling) is overtaken at low 
frequencies by the |cc;n|^~^0^ term associated with the 
long-range interactions of non-spin-flip type. 

To expand on this point, the short-range interaction 
term in the discrete Ising chain [Eq. (|4|)] that arises from 
mapping the transverse Kondo coupling can be expressed 



1 /'^^o 
log( — 



L 

E 



(12) 
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In the continuum limit, the corresponding kinetic energy 
written in terms of {dS^/dr)'^ clearly requires regular- 
ization through a finite value of tq. The long-range in- 
teraction in Eq. (|4|) must also be regularized, either by 
reinstating 00 c in Xo^(^) entering Eq. ([8]) or through 
enforcement of \r — r'\ > Tc in Eq. ([9]). Finally, tq 
serves to regularize the dynamically generated Kondo 
scale, Tk ~ (I/tq) exp(— l/Fro), below which the quan- 
tum scaling occurs. This scale vanishes when tq (or Tc) 
vanishes, making the quantum critical behavior inacces- 
sible. These observations suggest that the continuum 
limit, necessarily taken in the quantum-to-classical map- 
ping, fails to capture the the topological (vortex-like) ef- 
fect encoded in the Kondo spin-flips. What results is a 
change in the temperature dependence of the spin self- 
energy, qualitatively similar to that arising in the clas- 
sical model from changing the maximum winding from 
M = 1 to M = 00. 

A concurrent study [lO] addresses the mapping of the 
dissipative spin-boson model using a Monte Carlo algo- 
rithm that explicitly takes the limit tq ^ of a classical 
Ising chain [2l| . Ref. focuses on the finite-size scaling 
exponents of the static magnetic properties, rather than 
the scaling of the dynamical susceptibility. Its results 
for the mapped classical action are compatible with, and 
complementary to, our own. Ref. [13] also speculates that 
a marginal coupling introduced by bosonic truncation in- 
validates the NRG results ^ for the spin-boson model. 
However, our NRG studies for the BFKM, carried out 
over an extended range of bosonic truncation parameters 
and covering more than 20 decades of temperature, show 
evidence, neither for evolution in the exponent x away 
from its interacting value of 1 — e towards the mean-field 
X = ^, nor for the line of critical points that would be 
expected in the presence of a marginal operator. We 
also note that the NRG conclusion for the BFKM that 
X = 1 — e = y cannot be vitiated by a dangerously irrel- 
evant coupling. The NRG result corresponds to a spin 
self-energy [see Eq. ([11])] M{T) = Xo^ - 1/x ^ 
In the regime ^ < e < 1 of interest, this temperature 
dependence dominates any T^/^ term potentially gener- 
ated by a dangerously irrelevant coupling, in a manner 
reminiscent of what happens in the Monte Carlo calcu- 
lations with no winding (M = 1). These considerations 
all point to the quantum-to-classical breakdown being a 
real phenomenon. 

In summary, we have addressed the quantum-to- 
classical mapping of the Ising Bose-Fermi Kondo impu- 
rity problem. The finite-size scaling for the spin sus- 
ceptibility of the mapped classical chain demonstrates 
an intriguing dependence on the winding of the long- 
ranged interaction under periodic boundary conditions. 
Only for infinite winding does one recover the expected 
mean- field behavior. The contrast between these scaling 
properties and those of the Bose-Fermi Kondo model sug- 
gests a breakdown of the quantum-to-classical mapping 



for Ising-anisotropic quantum dissipative systems arising 
from the manner in which Kondo spin flips are treated 
in the continuum limit that is taken in such a mapping. 
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